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模糊 序列 的 GM(1,1) 建 模 * 


陈 定 元 ， 钟 金 标 
(安庆 师范 学 院 数学 系 ， 安 徽 安庆 246011) 

摘 E: 本 文 在 不 改变 GM(1,1) 模型 建 模 机 理 的 前 提 条 件 下 ， 运 用 模糊 回归 理论 对 GM(1,1) 模型 进行 优 
化 ， 通 过 定义 三 角 模糊 数 的 左 、 中 、 右 距离 ， 将 模糊 序列 的 GM(1,1) 模型 的 求解 转化 为 线性 规 
划 模 型 的 求解 ， 并 通过 实例 加 以 论证 。 研 究 表明 模糊 序列 的 GM(1,1) 模型 保留 了 GM(1,1) 模型 
所 需 建 模 数据 少 和 具有 预测 功能 的 特点 ， 且 能 为 决策 者 提供 一 个 决策 区 间 。 
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1 引言 


灰色 建 模 作 为 灰色 系统 理论 的 重要 内 容 ， 其 在 社会 生活 的 诸多 方面 得 到 了 广泛 的 应 用 。 为 
了 提高 模型 的 可 靠 性 与 适用 性 ， 许 多 科技 工作 者 在 邓 聚 龙 先生 提出 的 模型 的 基础 上 提出 了 针 
对 性 较 强 的 模型 L717。 研究 表明 ， 原 始 序列 数据 的 模式 及 光滑 性 是 影响 模型 精度 的 两 个 主要 因 
素 ， 因 而 形成 了 两 种 提高 GM(1,1) 模型 精度 的 方法 : 对 传统 的 GM(1,1) 模型 进行 修正 使 之 适应 
原始 序列 数据 的 模式 ， 对 原始 序列 数据 进行 变换 以 改善 其 光滑 性 。 现 有 研究 成 果 很 少 涉及 模 
糊 序 列 的 GM(1,1) 建 模 。 只 有 文献 [8] 讨论 了 对 称 三 角 模 糊 数 的 GM(1,1) 建 模 问题 ， 但 求解 方 
法 繁琐 ， 使 用 范围 有 限 。 本 文 将 在 文献 [8] 的 基础 上 利用 模糊 回归 理论 对 GM(1,1) 模型 建行 优 
化 ， 实 现 模 糊 序列 的 GM(1,1) 建 模 ， 并 探讨 相对 简洁 的 求解 方法 。 


2 GM(1,1) 模型 建 模 概述 
GM(1,1) 模型 主要 用 于 单个 时 间 序 列 的 动态 预测 ， 有 具体 过 程 如 下 ， 设 有 非 负 原始 时 间 序 列 
2 = {£0 (1), 192), ,z(o)(n)}. 
1) 对 原始 时 间 序 列 作 一 次 累加 生成 序列 zt = 4GCOz(0) 
z0 = {2 (1), 2(2),--- 2 (n)}, 
其 中 
(4) = > O(m), i=1,2,...,n. 


2) GM(1,1) 模型 的 灰 大 微分 方程 为 
ZO (Kk) +az(k) =b, (1) 
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其 中 i 
ZO (k) = 5 (eh) +29(k-1)), k=2,3,---,n, 


称 为 白化 背景 值 ，a 称 为 发 展 系 数 ，5 称 为 灰 作 用 量 。 其 对 应 的 白化 微分 方程 为 


dx) 
Q) 一 
ai + ax i (2) 


根据 最 小 二 乘法 原理 ， 式 (1) 的 参数 可 由 下 式 求 出 


| | = (BT B)-1BTYy, (3) 
—2z(2) 1 x!) (2) 
—z(3) 1 a) (3) 
_ Mn) l rO (n) 


对 应 于 式 (3) 求 得 参数 &,b， 则 式 (2) EMA 2M (1) = 101) 下 的 解 为 
4 (¢) = (zO0D - E eran 第 $. (4) 


3)” 作 一 次 累 减 还 原 2(0)(t) = TAGOR 的 便 可 得 到 原始 序列 的 拟 合 函数 


“is 2) (1) = 2 (1), t=1, 
E(t) = . ae 
et) -20(¢-1)=(1- eê) (2 (1) 一 ae t>1, 


44t=1,2,3,---,nit, (4) 就 是 式 (2) 的 离散 形式 解 
= PE Â\ alk- 
Ok) = (1 — e*) (2a) ~ z Je, (5) 


2 (k), k= 1,2,… ,n， 是 原始 序列 数据 的 拟 合 值 ，2(0) (k), kk > n， 是 模型 的 预测 值 。 
从 (5) 式 可 以 看 出 ，GM(1,1) 模型 的 模拟 与 预测 精度 取决 于 最 小 乘法 求 取 的 参数 6, 6 以 及 初 
值 z(o)(1)， 且 初 值 z(0(1) 对 生成 序列 zt0 和 原始 序列 的 模拟 预测 值 2(0) 均 具 有 指数 效应 四。 为 
了 消除 初 值 z(o)(1) 对 原始 序列 预测 值 的 指数 效应 ， 结 合 (4) 式 和 (5) 式 ， 可 以 将 预测 序列 的 表 

现形 式 优化 为 
2 (k +1) = (1 -= ê) (2 - = Jet App dia): (6) 


式 (6) 其 实 是 一 个 线性 回归 模型 。 根 据 模糊 线性 回归 理论 ，4o, Ai 不 仅 可 以 取 实 数 ， 而 且 可 以 
取 模 糊 数 。 当 Ao, Ar 取 模 糊 数 时 ， 式 (6) 表达 为 


zk+1)= Ap + Az (k). 


此 时 不 管 原始 序列 是 实数 序列 ， 还 是 模糊 序列 ， 输 出 序列 均 为 一 模糊 (区间 ) 数 序列 ， 从 而 
为 决策 者 提供 一 个 决策 区 间 ， 在 实际 中 具有 广泛 地 应 用 。 
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3 ”原始 时 间 序 列 为 模糊 数 的 GM(1,1) BR 
为 了 简化 模型 的 研究 ， 以 下 模糊 数 均 取 为 三 角 模 糊 数 。 当 原始 输入 序列 为 模糊 序列 ， 即 


z0) 一 {@ 1), ZO (2), ---, ZO (n)}, 


其 中 
FOG) = (ZO A)", (4), FO)®), 一 1 2 ,7 


原始 序列 通过 一 次 累加 得 到 模糊 生成 序列 为 ZO) = AGOZ 
20) = {#(1), #0) (2), +++, #™ (n)}, 


其 中 


EO (i) = (2 (a), EDS, z0 (i)®), #0) (3) 3 ZO (m), 


z0) (i)° = > ZO (m)®, EY G)FR = 3 g0) (m)®, i= 1, 2, en, 
m=1 


则 基于 回归 理论 的 GM(1,1) 优化 模型 为 
XO(k+1)= Ap + Až (k), (7) 


AF Ay = (Ab, AS, AÈ), Ar = (AZ, AF, AB), He AL, AS, 4AR 分 别 为 4;, i = 0,10 FR 
值 、 中 心 值 、 上 界 值 ， 则 它们 的 隶属 度 函 数 为 


a 
Z 一 45 
HAi(z) 1 TRAT AF <z < AR, 
0, 其 他 . 


运用 Zade 的 扩张 原理 ， 可 以 得 到 模糊 输出 的 隶属 度 函数 为 
z—(AS+ASa 1) (k)) c c 
l- ag aptat apog (48 + Aja (k)) 
<z < (Af + Agr) (k)) + (AR + Ağ) (k)), 
$ 2 (Ap +Aga™ (k))— 
MRO) =) 1 a apta ataa (48 + Af (k)) 
— (AE + Ar (k)) < x < (AS + Afx™ (k)), 
0, 其 他 . 


根据 模糊 回归 理论 ， e A1， 使 模糊 输出 的 隶属 度 大 
于 等 于 及 ， 且 左右 扩展 值 尽 可 能 小 。 


HXO RL) Z Ms (8) 
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其 中 及 (0 < h < 1) 由 用 户 确 定 ， 用 来 衡量 模型 对 样本 的 拟 合 程度 。 文 献 9] 建议 hh 值 选择 应 根 
据 原 始 序列 (PEA) 的 个 数 确 定 。 当 原始 序列 的 容量 充分 大 时 ， 可 以 选择 h = 0， 随 着 样本 容 
量 的 减 小 ，h 的 值 随 之 增 大 ， 但 中 心 值 不 变 。 文 献 [10] 推荐 的 取 值 范围 为 0.5 - 0.7. MA 
究 成 果 多 是 将 (8) 式 转化 为 线性 规划 问题 求解 ， 其 求解 方法 见 文献 [8]。 在 此 ， 针 对 模糊 序列 
的 GM(1,1) 模 型， 本 文 提出 一 套 相 对 简洁 的 计算 方法 。 

RRA XO (k +1) = o + Are (k) 的 求解 ， 根 据 文献 [11] 的 研究 思想 ， 引 入 模糊 数 左 、 
中 、 右 距离 的 定义 。 

定义 ” 设 三 角 模糊 数 A= (A>, AC, AR), Š = (Br, BC, BR)， 定 义 三 角 模 糊 数 距离 : 

1) 左 距 离 : d,(A, Š) = (AY — BŁ}; 

2) ÞES: dc(4, B) = (AC — ee 

3) AES: dp(A,B) = (AP — BR). 

实数 可 以 看 成 是 下 界 值 、 上 界 值 、 中 心 值 相等 的 特殊 模糊 数 。 根 据 三 角 模 糊 数 左 、 中 、 右 
距离 的 定义 和 性 质 以 及 回归 模型 的 求解 原理 ， 求 解 (7) 式 中 的 模糊 参数 Ao, A 就 是 其 在 模糊 数 
左 、 中 、 右 距离 的 意义 下 数据 


(XO F, OAE), (XO, OG), (KOR, ZO"), i=1,2,...,n—1 


的 误差 平方 和 最 小 。 根 据 左 右 扩展 值 尽 可 能 小 的 原则 ， 即 在 模糊 左 、 中 、 右 距离 的 意义 下 模糊 
拟 值 与 原始 序列 的 左 、 中 、 右 距离 最 小 。 所 以 求解 模型 羡 0(k 十 1) = Ap + Ape (k) 就 转化 为 


n-1 
min dz (4, Ap + Ait) = min $` (2 (+1) — (AE + Aba (i)"))’, 
i=l 


n—-1 
min do (€, Ap + Aya) = min 》 (2 (i + 1)° — (AF + AS 2 (4)°))?, 


i=1 


n-1 
min dp(#), Ap + Aix) = min X` (sO (i + 1)? — (AR + AP x) (4)®))?. 


i=l 


根据 最 小 二 乘法 ， 令 


n-1 
Bae = FeO + 1) — (Ag + Alo 4) =0, 
240 i=l 
adr “去 (0) L L L, (1) (7) L) (1) /NL 
FE = DOG +1)! = (Ak + Aba aH") =0 
1 i=1 
adc © (0) C L Cl) fC 
aac = D_ (eit DY — (Ag + Ape’ ()")) = 
OAG = 
ade n—i 


= y (zol 419° — (AL + AFV CWC) =0, 


i=l 


aAc 


400 


L 


adr n—-1 


OAR 


i=1 


求解 这 一 线性 规划 问题 得 


Ag = (AP, AŞ, AB), A, = (AP, AE, AÑ), 
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= > (zGTIR — (AF + APe™()*)) = 0, 


n—i 
QR Soot DR — (AB + AP ZO (iR) (@)*) = 0, 
i=l 


K(k +1) = Ay + Aye (k) = {A}, AÇ, AR} + {A}, AS, APE (Kk) 


= {XO (k +1), XOR41)°, XO (k +17}, 


其 中 XO)(k +1), XO), XO EDAX O FRA Poti EA 
值 。 根 据 (8) 式 ， 对 于 任意 的 he [0,1], & 


1O (k +1)? = AG — (1—h) AX + [AF — (1 — h) AE] 2 (k), 


2 (k +1) = AS + AS (Kk), 


t%(k+1)F = AY + (1— AAP + [AP +0- 


h)AP]2 (k), 


[2O@+1)], © [XOG@+1)],, Vee (0,1, #=1,2,---,n-1, 


与 文献 [12] 的 研究 结果 相 一 致 。 
算 例 文献 [8] 用 优化 GM(1,1) 模型 讨论 了 台湾 地 区 2001-2004 年 LCD TV 模型 销售 量 的 
GM(1,1) 建 模 问题 ， 现 用 其 数据 检验 本 文 的 关于 模糊 序列 GM(1,1) 模型 及 算法 的 可 靠 性 ， 并 加 


以 分 析 对 比 。 


利用 本 文 所 建 的 模糊 序列 的 GM(1,1) 模型 并 结合 上 述 算法 , 得 到 台湾 地 区 2001-2005 年 LCD 
TV 销售 量 的 模糊 模拟 与 预测 结果 ， 见 表 1。 将 本 文 模型 及 算法 结果 ， 与 收集 到 的 原始 序列 数 
HB, 文献 [8] 的 研究 结果 进行 对 比 ， 如 图 1 所 示 。 


表 1: ”本文 的 模糊 序列 的 GM(1,1) 模型 及 其 算法 的 可 行 性 分 析 


预测 值 

X (2001) 
X (2002) 
X © (2003) 
X © (2004) 
X (2005) 


统计 数据 

(73,81,89) 

(138, 150, 162) 

(377, 393, 409) 

(850, 870, 890) 
(2010.2, 2011.4, 2014.8) 


本 文 模型 及 算法 结果 
(73,81,89) 

(159.9, 169.9, 180.2) 
(347.1, 365.4, 383.6) 
(858.0, 877.6, 897.2) 
(2029.62, 2073.72) 


文献 [8] 模型 及 算法 结果 (h = 0) 
(73,81,89) 
(150.33, 194.43) 


(349.83, 393.93) 
(872.52, 916.62) 


通过 表 1 和 图 1 可 以 看 出 ， 本 文 所 建 模型 及 所 提 算 法 的 结果 与 收集 数据 和 文献 [8] 的 研究 结 
果 较 吻合 ， 因 此 是 可 行 的 。 
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4 


Se 收集 数据 下 界 值 
一 使 一 收集 数据 中 心 什 
一 。 一 收集 数据 上 界 值 
一 十 一 本 文 模型 的 下 界 值 
Se 本 文 模型 的 中 心 值 
E 本 文 模型 的 上 界 值 
一 雄一 文献 [8] 的 下 界 值 
—O— 文献 [8] 的 中 心 值 


0 
2001 2002 2003 2004 2005 


图 1: 数据 对 比 图 


结束 语 


本 文 根 据 GM(1,1) 模型 与 模糊 回归 模型 建 模 的 数学 机 理 ， 对 GM(1,1) 模型 进行 优化 ， 实 现 


了 模糊 序列 的 GM(1,1) 建 模 ; 根据 定义 的 三 角 得 模糊 数 的 左 中 右 距离 公式 ， 得 到 了 模糊 序列 
的 GM(1,1) 模型 的 一 套 相 对 简洁 的 算法 ， 通 过 实例 论证 了 模型 及 其 算法 具有 较 好 的 可 靠 性 ， 同 


时 ， 


模型 保留 了 传统 GM(1,1) 的 所 需 建 模 数据 少 和 具有 预测 功能 的 优点 ， 又 可 以 通过 模糊 输出 


为 决策 者 提供 一 个 决策 区 间 ， 具 有 一 定 的 实际 意义 。 


参考 文献 : 


[1] 


[2 


4 


邓 聚 龙 . 灰 理 论 基础 [M]. 武汉 :华中 科技 大 学 出 版 社 ，2002 

Deng J L. Grey Theory[M]. Wuhan: Huazhong University of Science and Technology Press, 2002 

WH, Xl AWE, 一 类 离散 灰色 模型 及 其 预测 效果 研究 [J]. 系统 工程 学 报 ，2006, 21(5): 520-523 

Xie N M, Liu S F. Research on discrete grey model and its predictive result[J]. Journal of Systems Engi- 
neering, 2006, 21(5): 520-523 

WAH, PB. GM(1,1) 模型 的 改进 方法 与 应 用 [可 . 系统 工程 与 电子 技术 ，2002, 24(4): 60-63 

Xu X L, Luo J. Improvement to GM(1,1) and its application[J]. Systems Engineering and Electronics, 2002, 
24(4): 60-63 


] 罗 党 ， 刘 思 峰 . 灰色 模型 GM(1,1) 优化 [可 . 中 国 工程 科学 ，2003, 5(8): 50-53 


Luo D, Liu S F. The optimization of grey model GM(1,1){J]. Chinese Engineering Science, 2003, 5(8): 50-53 
BA. 优化 灰 导 数 白 化 值 的 无 偏 灰 色 GM(1,1) 模型 [J]. 数学 的 实践 与 认识 ，2003, 33(3): 13-16 

Mu Y. An unbiased GM(1,1) model with optimum grey derivative’s whitening values[J]. Mathematics in 
Practice and Theory, 2003, 33(3): 13-16 

Wu C M, Zhang T, Gao X P. On discount weight grey modeling of GM(1,1){J]. Journal of China Universities 
of Posts and Telecommunications, 1999, 6(2): 78-82 

张 辉 ， 胡 适 耕 . GM(1,1) BAA AR A [J]. 系统 工程 理论 方法 应 用 ，2001, (1): 71-73 

Zhang H, Hu S G. Accurate solutioon for GM(1,1) model[J]. Systems Engineering-Theory, Methodology & 
Application, 2001, (1): 71-73 


402 I E 数 学 学 R 第 27 卷 


[8] Tsaur R C. Forecasting analysis by using fuzzy grey regression model for solving limited time series data[J]. 


Soft Comput, 2008, 12: 1105-1113 
[9] Tanaka H, Watada J. Possibilistic linear systems and their application to the linear regression model[J]. 


Fuzzy Sets and Systems, 1998, 27(3): 275-289 
[10] Bardossy A, Bogardi I. Fuzzy regression in hydrology[J]. Water Resources Resources Research, 1990, 26(7): 


1497-1508 

[11] WA, REB. 含 三 角 模 糊 参 数 的 线性 回归 模型 的 新 解法 [可 . 汕头 大 学 学 报 ，2004, 19(4): 1-9 
Cao B Y, Zhu Z X. A new method to the linear regression model with triangular fuzzy parameters[J]. 
Journal of Shantou University Natural Science Edition, 2004, 19(4): 1-9 

[12] Tanaka H, Uejima S, Asai K. Fuzzy linear regression model{J]. IEEE Transaction on SMC, 1982, 12(6): 


903-907 


GM(1,1) Modeling on the Fuzzy Sequence 
CHEN Ding-yuan, ZHONG Jin-biao 


(Department of Mathematics, Anqing Teachers College, Anging, Anhui 246011) 


Abstract: Without changing the principle of the GM(1,1) modeling, this paper adopts the fuzzy 
regression model to optimize the GM(1,1) model. By defining the left distance, the middle distance 
and the right distance of a triangular fuzzy number, the GM(1,1) modeling on the fuzzy sequence is 
resolved by changing it into a linear programming problem and demonstrated through an example. 
Research shows that the GM(1,1) modeling on fuzzy sequence not only keeps the character of needing 
few data and prediction, but also can provide a decision interval to the decision maker. 

Keywords: grey model GM(1,1); fuzzy regression model; fuzzy sequence; triangular fuzzy number; 
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